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^ ! Abstract 

^ _ The anomalous time depending blueshift, the so-called "Pioneer anomaly" , 

' that was detected in the radio-metric data from Pioneer 10/11, Ulysses and 

CN . Galileo spacecraft may not result from a real change of velocity. Rather, 

I the Pioneer anomaly may be understood within the framework of general 

■ relativity as a time depending gravitational frequency shift accounting for the 

. time dependence of the density of the dark energy when the latter is identified 

with quintessence. Thus, instead of being in conflict with Einstein equivalence 
principle, the main Pioneer anomaly appears merely as a new validation of 
general relativity in the weak field and low velocity limit. 
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1 Introduction 

Since 1998, Anderson et al. have continuously reported an anomalous Doppler shift 
derived from a twenty years study of radio-metric data from Pioneer 10/11, Ulysses 
and Galileo spacecraft [1]. The observed effect mimics a constant acceleration acting 
on the spacecraft with magnitude ap = (8.74 ±1.33) x 10"^ cm s~^, directed toward 
the sun to within the accuracy of the Pioneers' antennas and a steady frequency 
drift ~ 6 X 10-9 Hz/s which equates to a "clock acceleration" ^^^^ = . An 
independant analysis of radio Doppler tracking data from the Pioneer 10 spacecraft 
for the time period 1987 - 1994 confirms the previous observations [2]. Besides, it 
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has been noted that the magnitude of op compares nicely to the Modified Newto- 
nian Dynamics (MOND) acceleration constant oq ~ cHq [3], where Hq is the Hubble 
parameter at present cosmological epoch. Actually, the reported anomaly cannot 
really be related to MOND. Indeed, the gravitational pulling of the sun up to 100 
AU still yields an acceleration greater than qq by at least three orders of magni- 
tude, equating oq only at about 3000 AU. Hence, Newtonian dynamics up to GR 
corrections should apply to spacecraft. The same argument was put foreward by 
M. Milgrom [4] to reject the claim that MOND fails in the laboratory [5]. In a 
further study, Anderson et al. found that the small difference of magnitude of ap 
for Pioneer 10 and Pioneer 11 is related to their difference of spin-rate history [6]. 
Removing the spin-rate change contribution, yields an apparent anomalous acceler- 
ation ap = (7.84 ± 0.01) X 10~^ cm of the same amount with a great accuracy 
during a long time interval for both Pioneer 10/11. Since this has found a con- 
ventional explanation [7], it points out in favor of an external origin for the main 
Pioneer anomaly at the expense of any possible internal cause (e. g., heat system- 
atics). As yet, this also excludes any significant drag force acting on the spacecraft 
(e. g., the mirror world explanation [8]) as an explanation of the anomaly. Indeed, 
the velocities of the Pioneer 10 and 11 relative to the sun differ by about 5% in 
the time interval of interest. Moreover, Whitmire and Matese [9] have carried out a 
statistical analysis by comparing the mean original energy obtained from the Comet 
Catalogue of Marsden and Williams [10], with and without corrections due to an 
hypothetical Pioneer anomalous acceleration, to that predicted by galactic tidal the- 
ory. They conclude that the implied higher binding energy is incompatible with the 
established evidence that the galactic tide is dominant in making Oort cloud comets 
observable. The whole work of these authors is based on the major assumption that 
the Pioneer anomaly reveals a real acceleration. As they emphasized, systematic 
outgassing would not hide the efi'ects of such anomalous acceleration if real. There- 



2 



fore, the alternative that the Pioneer anomaly does not result from a real change in 
velocity deserves to be investigated. Indeed, a direct interpretation of the observa- 
tional data from the spacecraft implies merely a time dependent frequency blueshift 
of the photons. However, any true Doppler shift would involve an accompanying 
acceleration. Hence, given our knowledge of photons frequency shift, the only other 
relativistic effect that can be confused, at the solar system scale, with a real Doppler 
shift is the gravitational frequency shift. In the weak field and low velocity limit, 
this should involve a time dependent gravitational potential instead of a spatial de- 
pendent one, in order to avoid an induced anomalous acceleration for test bodies. 
Of course, the origin of a time dependent gravitional potential would not be relevant 
if not justified on the physical ground. Actually, although unfamiliar, we show that 
a time dependent gravitational potential follows from our present knowledge of the 
matter-energy content of the universe. Indeed, quintessence is one possible form of 
dark energy that explains the recent discovery of the accelerating universe in term 
of a scalar field whose energy density is time dependent and dominant over that of 
the ordinary matter at present cosmological epoch. Now, a time dependent matter- 
energy density will naturally generate a time dependent potential. Given that the 
energy density of the quintessence is on average the same everywhere depending only 
on time, this should holds even at small scale, not only at the cosmological level. 
The idea of a time dependent gravitational potential have been revived recently by 
A. F. Ranada [11] and K. Trencevski [12]. However, both attempts introduce only 
phenomelogically the time dependent potential (e. g. by choosing its sign by con- 
venience and trying to match exactly ap with cHq). As a common feature, none 
of both approaches refer to a field equation to derive the time dependent potential. 
Above all, none rely intirely on GR. Besides, we note that the link between the time 
dependent potential and dark energy (whether this pertains to a true cosmological 
constant or quintessence is not expressed at all) is either put by hand [11] or just 
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invoked by words [12] by the authors. The aim of this paper is to demonstrate that 
the Pioneer anomaly can find a natural explanation within the framework of gen- 
eral relativity (GR) when the gravitational contribution of quintessence is properly 
accounted for. Then, the Pioneer anomaly interprets as a gravitational blueshift 
in excess induced locally by the energy density of quintessence which provides a 
time dependent contribution to the gravitational potential. Furthermore, because 
of its weakness, the Pioneer anomaly seen as a time dependent gravitational po- 
tential actually entails no contradiction with respect to the four classical tests of 
GR. Also, it is worth noticing that the analysis that is carried out in the following 
involves the Hubble constant only through the time dependence and non locality of 
the energy density of the quintessence. In any case, the matter content of gravita- 
tional bound systems like galaxies or the solar system arc almost decoupled from 
the general expansion of the universe whose effect is quite negligible on the stars 
or planetary orbits (see [13] and [14], section 16.1 H). Indeed, the exact mathemat- 
ical treatments of gravitational bound systems in an expanding universe that are 
based on GR solely (with or without a true cosmological constant) always lead to 
conclude to a negligible effect out of an ensemble of dust particles [15, 16, 17] or a 
star [18]. Thus, although the net effect is not exactly zero, it cannot account for the 
Pioneer anomaly in the framework of GR (see e. g., [19]) without reconsidering the 
matter-energy content of the universe. 

2 Time varying potential from quintessence 

Let us consider the dark energy as quintessence, that is a neutral scalar field cou- 
pled to ordinary matter only through its gravitational infiuence. We show that this 
involves a time varying contribution to the gravitational potential even for gravita- 
tionally bound systems. As one knows, unlike the ordinary matter (luminous or dark 
matter) whose mass density tends to decrease continuously driven by the expansion 

4 



of the universe or clumps locally under the gravitational pulling, the energy density 
of the dark energy remains constant on average everywhere at any given cosmic time 
(quintessence models) if not a universal constant (a property of spacetime itself) in 
the case of a true cosmological constant. For that reason, we will not assume a static 
metric for the solar system hereafter. Instead, in addition to the point like mass of 
the sun and the other heavenly bodies (whose effects as usual are seen as perturba- 
tions), we will also deal with the gravitational effect of the uniform mass density, 
Pq, of the dark energy which behaves like a cosmological "constant" A = SttGpq/c^. 
By choosing the coordinates such that the metric tensor be diagonal, the metric still 
expresses in the canonical form 

ds' = ( 1 + 2 ^) c^de - ( 1 + 2 ^ )-i dr' - ( ^ sin^ ), (1) 

but with the potential V now depending ton both spatial and time coordinates, 
where — ct, — ris the radius from the sun, x'^ — 9 is the polar angle and x^ — (f 
is the azimutal angle. The function V is derived from Einstein equations i?^,^ — 
Qfiu — ill the limit of the weak and slowly varying gravitational 

field. We proceed as usual by setting g^^, — rj^i, + h^,^, where all the /i^^i^'s are much 
less than unity. Thus, the Einstein equations reduce to Rqq — ^^{pmatter+PQ+^^) 
in the weak field and low velocity limit. Expanding Qqq and gu in the first order 
in V/c^ and Rqo in the first order of the derivatives of the h^^^s (see e. g., [14], 
section 15.1 and [20]) but without discarding as usual the time derivatives, one gets 
Roo = I ( V,oa - ^oo,/3a + ^oa,/30 - ^/3a,oo ) = ^{^^ + AF). Heuce, the 
equation to solve reads 

+ AV^AnG (^pmatterir) + pQ + 3^y (2) 

where the equation of state of the quintessence is given by Pq = wq pqc^ with 
equation of state parameter wq ~ — 1. Usually, one solves (2) by assuming a static 
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potential. This yields (up to an integration constant) 



V = Vn 




(3) 



where the Newtonian potential Vn = V/v(r) is a solution of Gauss equation 



and the quadratic extrapotential term is related to the dark energy density. How- 
ever, because of the fine-tuning that would be required both for the cancellation 
of the vacuum energy (quantum field theory problem) and the coincidence prob- 
lem (why the density of the dark energy and that of the ordinary matter are of 
the same order just today, whereas both rates of evolution are quite different), 
quintessence models are prevailing in the literature at the expense of a true cos- 
mological constant [21]. A quintessence model involves a dynamical cosmological 
constant driven by an almost homogeneous scalar field (spatial fluctuations cancel 
out on average) [22], Q, with a suitable self-interaction potential (most often chosen 
as an inverse power law of Q). This means that in some sens, though not exactly, A 
is deflned as a function of time in such models. So, the previous solution (3) which 
is purely static does not hold any more. Thus, while in the case of a true positive 
cosmological constant the vacuum solution yields a de Sitter spacetime, spacetime is 
rather asymptotically Minkowskian in the case of quintessence whose energy density 
decreases monotonously toward zero with respect to time. Since pq — pgit) is non 
zero everywhere, Birkhoff's theorem cannot rigorously be invoked to infer a static 
spacetime. Let us emphasize that the generalization of Birkhoff's theorem to the 
case with a A-term (cosmological constant) is restricted to the case when the latter 
is indeed a true universal constant [23]. Besides, it is clear that making a gauge 
transformation in order to cancel out any time dependent potential term will not be 
hepful, since this is associated with a change of reference frame. Now, we need to 



AVat = 47rGpmatter 



(4) 
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study the motion of spacecraft in the same reference frame as for the planets. So, 
let us look for a solution of the form V — VN{r) + VQ{t). Equation (2) then splits 
into equation (4) and the following 

^^=-S.Gp^(ty. (5) 

At present cosmological epoch, pQ{t) is a slowly varying function as compared to 
the time scale that concerns the Pioneer 10/11 spacecraft. Therefore, the solution 
of equation (5) reads in the first approximation 

VQ{t) = - ^ ilQHyit - to)'' - mQHlc\h - m - to) + VQ{to), (6) 

where ti is the cosmic time at which VQ{t) passes by its maximum value and 
Qq — A{to)c^/SHQ denotes the density parameter associated with the quintessence 
at present cosmological epoch. Let us point out that one passes from the non- 
Newtonian extra potential of relation (3) to relation (6) by making the substitution 
r — > c(i — ti) in (3) then multiplying the whole result by 3 (because of the isotropy 
of the 3-space) . As regards the phenomenology, it seems reasonable to assume that 
the linear term of eq(6) is of the order of the quadratic one early in the past but fi- 
nally dominates over it at present epoch. These conditions are satisfied for ti ~ to/2 
and imply VQ{to) ~ Vq{0). Present cosmological observations yield fig ~ 0.7 and 
Hoto ~ O.QGjloloe for a flat universe [24]. Hence, relation (6) yields 

Ym^-H^^it-to)'-Ho{t-to) + (7) 
c c 

According to recent observations [25], the cosmic jerk when the expanding universe 

made the transition from deceleration to acceleration corresponds to a redshift Zj — 

0.46 ± 0.13 or equivalently a fractional look-back time {to — tj)/to ^ 0.4. So, the 

choice of parameter made above is consistent with the natural expectation that the 

potential energy Uq = rnVg of a test particle of mass m decreases monotonously 

past tj. 



3 Interpretation of the observational data 
3.1 The main Pioneer anomaly 

Since g'oo = 1 + 2 ( ^^^^"^ ) , the resulting gravitational frequency shift of photons is 
derived from the following relation 



— ~r ^ordinary ~r ^ anomalous ^ ) 



T^A \ 9oo{B) 

where Zordinary — ^ ^"^^ represents the familiar gravitational frequency shift 
whereas Zanomaious = ^'^^^^^ is the time varying gravitational frequency shift 
related to the quintessential dark energy. For a small time of flight, At = Ib — tA, 
of the photons of the communication signals between the spacecraft and the earth 
(At negligible as compared to the age of the universe), one gets in the first order 
approximation on account of relation (7) 

1 dVQ , tyl + t^ s . /„s 

^anomalous — ^ ~~dt~ 2 (y) 

Since ^ < 0, relation (9) imphes a systematic blueshift in excess. For the time 
interval of interest for the Pioneer 10/11 spacecraft, the quadratic term of relation 
(7) is negligibly small with respect to the linear one. On account of relation (9) 
above, this involves an anomalous blueshift in excess 

^ Op At 



^anomalous 



(10) 



c 

where ap Hqc in accordance with observations. 

The above analysis can be pushed further by taking into account the revolution 
and rotation of the earth. Then, one derives the annual and dayly modulations 
following the same reasoning as before. We argue that these modulations are fur- 
ther evidences for the dependence of the spacecraft anomalous frequency shift with 
respect to time, in particular the time of flight of the photons of the communication 
signals. 



3.2 The absence of a true anomalous acceleration 

Moreover, the extrapotential Vq will not effect the motion of the spacecraft or heav- 
enly bodies being time dependent only. Indeed, in the weak fields and low velocity 
limit, the geodesic equation simplifies to 

dV _ 1 

= -Wiv, (11) 
where r denotes the position vector of the test body. 

3.3 The steady frequency drift 

According to Anderson et al. (see [6], section V - B and C, p. 19), the "clock 
acceleration" at — d{Au/u)/dt cannot be attributed to a systematic drift in the 
atomic clocks of the DSN (Deep Space Network, the network of ground stations 
that are employed to track interplanetary spacecraft). We show that it follows from 
the path equation of photons dP'^ — — F^^ P" dx^ for the time component — 0, 
where po = ^ and P — define the 4-momentum P^ of a photon of frequency u 
propagating in the direction of the unit vector u. Since Fqq ~ ^ ~ — ^ in weak 
field and low velocity limit, one finds after removing the gravitational frequency shift 
due to the Newtonian potential solely 

1 , du . du . . ap 

as observed. Let us notice that relation (12) can be derived in the usual simple 
manner from the conservation of energy for the photon to which one assigns a 
gravitational mass m — ^ in accordance with the EP and a kinetic energy Ec — hv. 
Further, the other first order correction implies by the potential Vq comes from the 
radial component = 1. This brings an additional contribution to the deflection 
angle of the photons near the sun of the order A(5 ~ at for a time of flight Ai. 
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Clearly, AS is by far less than the classical result S = 1.75" by more than 7 orders 
of magnitude. 

4 Discussion 

The study that we have carried out in this paper does not pretend to bring the last 
word on the so-called Pioneer anomaly. We have just tried to clarify the subject as 
to the consistency of such an effect with general relativity. Although, the results 
seem encouraging in this respect, we still have at least two problems to face. Namely, 
whereas the slope of the observed anomalous frequency shift versus time appears 
almost constant between 20 AU and 80 AU (see e. g., FIG. 1 of reference [1] for 
Pioneer 10 in the time period 1987-1995), it drops down to zero below 10 AU (see 
FIG. 6 or FIG. 7 of reference [6]). Nevertheless, Ulysses and Galileo which where 
moving below or about 5.4 AU (orbital radius of Jupiter) also have shown an almost 
apparent constant acceleration in their radio Doppler and ranging data but with 
a high correlation with solar pressure (respectively 0.888 and 0.99). These issues 
have never been addressed up to now by any theoretical study. Of course, the latter 
observations cannot fit within the framework of GR as a result of quintessence since 
the EP would be strongly violated. As a way out, we conjecture that mass loss 
toward the sun may potentially explain, at least in part if not all, the vanishing 
of ap observed on both Pioneer 10/11 spacecraft below 10 AU. Indeed, since the 
temperature within the Pioneer 10/11 is almost stabilized to about 300 K while 
decreasing with respect to the heliocentric radius, outgassing out of the spacecraft 
may occur (note the large error bars on ap below 10 AU for all four spacecraft). 
Also given that the mass of the Pioneer has decreased by almost 7% since the date of 
launch, a mass loss of about 2 ppm is still possible during the following ten years after 
launch. Because of the gravitational pulling of the sun, we expert that the ejected 
mass will be directed toward the sun. For both Pioneer 10/11, one just needs an 
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almost steady rate of mass loss directed toward the sun m ~ 0.5 gm per year during 
the ten years following the launch of the spacecraft but quickly decreasing past this 
time interval. As for the issue of the high correlation with the solar pressure in 
the case of Gahleo and Ulysses, given the proximity of the sun, it is hkely that the 
capture of dust particles of the solar wind may provide an answer. Let us recall 
that the Pioneer 10/11, Galileo and Ulysses missions were also dedicated to collect 
dust particles from the solar wind, all the four spacecraft being endow with dust 
detectors on-board [27]. A more detailed study including precise estimates of the 
aforementioned effects will be presented elsewhere. 

5 Conclusion 

The study of the Pioneer 10/11 spin-rate change histories points out in favor of the 
possibility that the cause of the main Pioneer anomaly is external to the spacecraft. 
However, interpreting the main Pioneer anomaly as a true acceleration leads to some 
difficulties as regards its effect on Long-Period Comet Orbits and the Oort cloud. 
Especially, assuming a resisting medium as cause of the anomaly (with the drag 
acceleration proportional to the orbital speed or its square) not only necessitates a 
hole inside 10 or 20 AU to avoid large effects on short period comets such as Halley 
and Encke [26] (note also how huge such an anomalous acceleration would be relative 
to the earth-moon system, if real about the orbital radius of the earth) but above all 
also makes difficult to understand the analogous anomalous accelerations observed 
on Ulysses and Galileo data (orbital radius of both spacecraft less than 5.4 AU). The 
alternative view that we have adopted then consists to consider that the Pioneer 
anomaly may not be the result of a new force term. Indeed, the main Pioneer 
anomaly may be well understood within the framework of GR as a time depending 
gravitational frequency shift (with respect to the reference frame used by JPL for 
Doppler tracking of spacecraft) . The latter is derived from an extra potential whose 
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source is the time depending energy density of quintessence like dark energy Thus, 
the Pioneer anomaly may accomodate the EP as expressed by the relation at — ap/c 
between the clock "acceleration", a*, and the anomalous "acceleration", ap. Finally, 
the Pioneer anomaly seems to be rather a new validation of GR in the weak field 
and low velocity limit. Moreover, the interpretation of the Pioneer anomaly in the 
framework of GR favors clearly the case of quintessential dark energy at the expense 
of a true cosmological constant. This exhibits the potentiality of quintessence to 
help GR gives rise to some Machian behavior. The features emphasized above 
should make clear the very difference between our new approach and any other one 
published as yet on the subject. 
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